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Abstract 

The generating functional for Green functions of quark currents is given in closed form 
to next-to-leading order in the low-energy expansion for chiral SU{3), including one-loop 
amplitudes with up to three meson propagators. Matrix elements and form factors for 
strong and nonleptonic weak processes with at most six external states can be extracted 
from this functional by performing three-dimensional flavour traces. To implement this 
procedure, a Mathematica® program is provided that evaluates amplitudes with at most 
six external mesons, photons (real or virtual) and virtual (semileptonic form factors). 
The program is illustrated with several examples that can be compared with existing 
calculations. 
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1 Introduction 



Chiral perturbation theory [1-4] is the effective field theory of the standard model at low 
energies. In the mesonic sector, the state of the art is next-to-next-to-leading order in the 
low-energy expansion for strong processes [5] and next-to-leading order for nonleptonic weak 
transitions. 

Although most one-loop amplitudes for processes of physical interest have already been 
calculated a compact closed form for such amplitudes would still be useful, both in principle 
and in practice. The generating functional for Green functions of quark currents is the relevant 
object for this purpose. The generating functional to next-to-leading order was already calcu- 
lated in the classic papers of Gasser and Leutwyler [2, 3] for the strong interactions, including 
electromagnetic and semileptonic weak form factors, and by Kambor, Missimer and Wyler [6] 
for nonleptonic weak processes. In both cases, the calculation was limited to processes where 
at most two meson propagators occur in one-loop amplitudes. 

More recently, the generating functional for chiral SU{2) was obtained in closed form for 
transitions involving up to three propagators in one-loop amplitudes: at most three pions in 
Ref. [7] and at most two pions and a photon in Ref. [8]. The purpose of this article is to 
calculate the analogous generating functional with up to three meson propagators for chiral 
5*^7(3) for both strong and nonleptonic weak transitions. The corresponding SU (3) calculation 
with dynamical photons is in preparation [9]. 

Transitions involving at most three meson propagators in one-loop amplitudes comprise 
almost all processes of physical interest at next-to-leading order. More precisely, the complete 
strong and nonleptonic weak amplitudes and form factors to that order can be extracted 
from the generating functional for transitions of at most 0(0^), i.e. with at most six external 
mesons, photons (counting as 0(0^)) and virtual W^. Although most of such processes have 
already been calculated to next-to-leading order the generating functional provides a closed 
expression for such amplitudes with the same conventions and notation. As experience has 
shown, such a presentation greatly facilitates comparison with previous work and eventual 
detection of misprints and other errors. In addition to reproducing previous results in a simple 
and straightforward manner, several nonleptonic transitions have been calculated for the first 
time [10] with the generating functional presented here. 

Although the extraction of amplitudes or form factors from the generating functional boils 
down to performing three-dimensional flavour traces such a procedure may turn out to be 
quite time consuming, especially for transitions with three-propagator loop contributions. To 
facilitate this work, we provide the Mathematica® [11] program Amp calculator, written by 
one of us (R.U.), that evaluates the necessary traces upon input of at most six external states 
(mesons, photons, virtual W^) with their corresponding momenta. Use of the program is 
straightforward and does not require detailed knowledge of Mathematica®. 

The generating functional is introduced in Sec. Eland its low-energy expansion is discussed. 
The generating functional of next-to-leading order is treated in Sec. 01 for the strong interac- 
tions. We extend the work of Gasser and Leutwyler [2,3] by calculating the explicit form of the 
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one-loop functional for up to three meson propagators and to at most 0(0^). The renormaliza- 
tion of the generating functional is briefly reviewed. In SeclUthe nonleptonic weak interactions 
are included. The additional terms in the generating functional linear in the weak coupling con- 
stants Gs, 6*27 are calculated under the same conditions as for the strong part. In Sec. El we 
discuss the Mathematica® program Ampcalculator for evaluating amplitudes and form factors 
for a given set of external states to next-to-leading order. The program is exemplified for the 
decays —>■ n^l^ui {Kjf^), K~ tt^ti^, —>■ n^'-f'-f and the results are compared with those 
in the literature. Our conclusions are summarized in Sec. El Details about the input for the 
program Ampcalculator, one-loop integrals and constituent functions and the strong and weak 
Lagrangians of 0{p'^) are collected in three appendices. 



2 Generating functional of Green functions 

The generating functional for Green functions of quark currents is defined in terms of the 
vacuum transition amplitude in the presence of external fields, 

^iZ[v,a,s,p] Q^^io in >^_^_^_^ ^ (2.1) 
associated with the Lagrangian [2, 3] 

^ = ^QCD + ^7^K + «/.75)g - q{s - ipi5)q , (2.2) 

where C^q-q is the QCD Lagrangian with three massless quarks and Vf^,af^,s,p are three- 
dimensional hermitian matrix fields. 

Green functions of quark currents are obtained by functional differentiation of Z[v, a, s,p] 
in the usual way. To implement explicit chiral symmetry breaking, the external scalar field s 
is set equal to the quark mass matrix A4 = diag(mu, m^, m^) at the end of the calculation. 
The external spin-1 fields generate electromagnetic and semileptonic form factors with the 
assignments 

r^' = v'' + = -eQAf" , (2.3) 

l^' = v^'-a^' = -eQA^' - ^J- {W^^T+ + h.c.) 

v2 sin Q^J^ 

for external photons (A^) and W bosons {W^\ Qw is the weak mixing angle, Q is the quark 
charge matrix and T+ contains the relevant elements of the Cabibbo-Kobayashi-Maskawa ma- 
trix: 

/ I \ / K,^ K,. 

Q = 



V 




3 " " \ " ud ^ us 

-| 1 , 



I . (2.4) 
V 



At the hadronic level, the generating functional may be calculated in terms of an effective 
Lagrangian of pseudoscalar mesons and external fields [4], 
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The mesonic effective chiral Lagrangian will be needed to next-to-leading order here, including 
both strong and nonleptonic weak parts: 



£efT {U, v,a,s,p) = C2 + C4 + . ■ . 
The strong chiral Lagrangian takes the well-known form [3] 

F2 



^2 



-{D^UD^'U^ + xU^ + X^U) , X = '^B{s + ip) 



(2.6) 



(2.7) 



in terms of the pion decay constant and the quark condensate in the chiral limit: 

= F[l + 0{mg)] = 92.4 MeV (2.8) 

{0\uu\0) = ~F^B[l + 0{mg)] . 

(. . .) denotes the 3-dimensional flavour trace. The strong chiral Lagrangian Cf of 0{p^) [3] is 
reproduced in App. [Oand the nonleptonic weak Lagrangian will be discussed in Sec. 0] 

The chiral expansion of Ccs induces a corresponding expansion for the generating functional: 



Z = Z2 + Zi + ... 
At lowest order, the functional Z2 is equal to the classical action 

Z2 [f , a, s, p] = / ci^x C2iU,v,a,s,p) 



(2.9) 



(2.10) 



where U = U[v, a, s, p] is to be understood as a functional of the external fields via the equation 
of motion (EOM) for C2 (given here for the strong interactions only) : 



(2.11) 



Mesonic amplitudes at lowest order can therefore be read off directly from the lowest-order 
Lagrangian ()2.7|) and the corresponding nonleptonic weak Lagrangian CY ()4.H) by us- 
ing an explicit parametrization of the matrix- valued meson field U{ip), e.g., the exponential 
parametrization 



\/2 ^ 



f/((p)=exp(a„^7F) , -^K^^ 



\ 



IT 



K- 



K+ 



\ 



2VS 

Vg J 



(2.12) 



The matrix field ()2.12|1 also defines our sign conventions for mesonic amplitudes generated by 
the Mathematica® program Ampcalculator to be discussed in Sec. El 
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3 Generating functional of O(p^) 



The generating functional at next-to-leading order consists of three parts: 

■^4 = Zf'^'^ + Z^~^ + 2^wzw • (3-1) 

The tree-level part Zf*^*^ is given by the action for the next-to-leading Lagrangian C4^{U, v, a, s, p) 
to be taken again at the classical solution U satisfying the EOM ()2.11|) . The Wess-Zumino- 
Witten functional ^wzw [12] accounts for the chiral anomaly. We will not reproduce its explicit 
form here but it will be implemented in the Mathematica® program Ampcalculator (cf. Sec. EI). 

The one-loop functional Z^^^ is calculated in the standard way. We repeat the main steps 
here in order to define our notation. The matrix field U{ip) is expanded around the classical 
solution U = in terms of a traceless, hermitian fluctuation matrix ^: 

U = u^ = u{l + i^-^e + ---)u. (3.2) 
Working with dimensional regularization in d dimensions, one obtains 



■yL=X r F' d^xCixW^X^Cix) 1 

e^^4 = jd^^[i\ e 2 J ^ ^ ^ ^ }^ ^ ) ^ ^ ^ _^^^a (3 3) 

The measure c?/i[^] is defined in such a way that Zf"^ vanishes when the external fields f , 
a and p are set to zero and s is set equal to the quark mass matrix. In the same limit, the 
differential operator 'D"'^ turns into the Klein-Gordon operator with the appropriate meson 
masses to leading order in the chiral expansion^. The one-loop functional can then be written 
in the form 

Zf=^ = -InDet:^ = -Trlu:^ . (3.4) 

Splitting the differential operator V into the Klein-Gordon operator Vq and a remainder 5, one 
gets 

Zf=^ = ^Trln^^ = ^Trln(l + i) = ^ Tr In (1 - 5A) . (3.5) 
2 Uq 2 2 

The diagonal matrix A is the inverse of —Vq and it contains the Feynman propagators of 
pseudoscalar mesons: 



[2iTpp'^ — + ^e 

For the strong interactions 6 is of the form [3] 

6s = {r^,d,} + r^^r, + a (3.7) 

^ Since we work in the isospin limit throughout this paper the differential operator Dq is diagonal. 
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with 

aab = \{[K.y^.][Xb.y^]) + \{{\aAb}(y) , (3.8) 

To include one-loop contributions with up to three meson propagators, the one-loop functional 
()3.5p must be expanded to third order in 5: 

Zf=^ = - J Tr (5A) - ^ Tr ((5A(5A) - J Tr (5A(5A(5A) + . . . . (3.9) 
2 4 6 

The functional ()3.9j) suffices for transitions up to 0(0^), i.e. with at most six external mesons 
where an external photon counts as 0{(fP'). We will therefore limit the further discussion to 
processes of at most 0(0^). Since the external fields and s couple to at least two pseudoscalar 
meson fields whereas and p couple to at least one field, the functional needs to be expanded 
to third order in and s and to sixth order in and p. Keeping in mind that the matrices 
and a are of 0(0^), the strong one- loop functional assumes the form 

Zf=i= Tr(a(x)A(0) + f'^(x)f^(a;)A(0)) + 



I 
4 



I 
6 



I 
2 



Tr ({f^(a;), a.JA(x - y){r (y), a,JA(y - x) + 

+2{f^(x), 9, JA(x - y)V-{y)VM^{y - ^) + 
+2{f^(a;), a, JA(x - y)a{y)^{y - x) + 

+2f'^(x)f^(x)A(x - y)a{y)/^{y - x) + a(a;)A(x - y)cT{y)/\{y - x)) + 

Tr ({f^(x), 9.JA(x - i/){^(y), 9,JA(y - z){V\z), 9,JA(z - x) + 

+3{f^(x), 9. JA(x - y){V^{y), 9,JA(y - ^)a(^)A(^ - x) + 
+3(t(x)A(x - y)a(y)A(y - ^){f^(^), 5,JA(z - x) + 
+o-(x)A(x - y)(f{y)A{y - z)a{z)A{z - x)) + . . . 

Tr (a(x)A(O) + f^(x)f^(x)A(0)) + 
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+t^{x)dyM^ - y)t'{y)d.,Hy - x) + 
-V'^{x)d^^dyM^ - 2/)^(2/)A(2/ -^) + 
-f'^(x)A(x - y)V''{y)d^^PyMy - ^) + 
+2f'^(x)9,.^ A(x - y)t''{y)t,{y)^{y - x) + 

-2f^(x)A(x - y)t''{y)t,{y)d^^A{y - + 
+2f''(x)9,^A(x - y)cr{y)/\{y - x) + 
-2f^(x)A(x - y)a{y)d^^A{y - x) + 
+2f''(x)f^(x)A(a; - y)a{y)/\{y - x) + 
+a{x)/\{x - y)a{y)/\{y - z)) + 

Tr {t^{x)d.,^A{x - y)t'{y)dyMy - z)t\z)d^A{^ - x) + 

-f'^(x)9,„ A(x - y)t"{y)d,My ' z)t\z)d^^/\{z - x) + 
+f^(x)9,„ A(x - y)r (y)A(y - z)t\z)d^^d-^,A{z -x) + 

+V>^{x)d^^dyM^ - y)t''{y)d-,A{y - z)V\z)/\{z -x) + 
+f^(x)A(x - y)t'^{y)dyA^A{y - z)t\z)d,^A{z - x) + 
-t^{x)A{x - y)t^{y)dy,A{y - z)t\z)d,^d,,A{z - x) + 
-T''ix)d,^dyM^ - y)r'^(y)A(2/ - z)T\z)d,,Aiz -x) + 
-t^^{x)d^^A{x - y)t''{y)dy^d,^A{y - z)t\z)A{z - x) + 
+3t^^{x)d,^A{x - y)t''{y)dyMy - zMz)A{z - x) + 

+3t^{x)dyMx - y)f"'iy)Hy - z)a{z)d,^A{z -x) + 

-3f^(x)A(x - y)T''{y)dy,A{y - z)cx{z)d,,^A{z - x) + 
-3f^(x)9,^9„ A(x - y)^(i/)A(y - z)aiz)Aiz - x) + 
+3a{x)A{x - y)(y{y)A{y - z)t^'{z)d,^A{z - x) + 
-3a{x)A{x - y)a{y)d,^A{y - z)t^'{z)A{z - x) + 

+a{x)A{x - y)a{y)A{y - z)a{z)A{z - x)) + . . . (3.10) 



Integration by parts was used in Eq. ()3.10|1 to shift the derivatives from the matrices containing 
the fields to the propagators. In terms of the functions A and Gi defined in App. |B1 the strong 
one-loop functional to 0(0^) can be written 




p 
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Pi, 



Figure 1: Definition of external momenta. 



P,Q 

+ fd'xd'y^^e-^^(^-y^ E {GrpQ{p)t,Qp{x)t^pQ{y) 
•' y'^'^i p,Q,R 

+G%q{p) t^Qp{x)apQ{y) 

+Gzpq{p) {aQp{x)apQ{y) + 2T pqr{x)TP Rp{x)apQ{y 
+ I d'^xd^yd'^zf^^-^^e'^-^-^-^h'^^^y-^^ 



X 



(27r)4 (27r) 

E (,<^4PQi?(Pa,P6) t^,Qp{x)f^pR{y)fxRQ{z) 

P,Q,R 

+G5PQR{Pa,Pb) T^^Qp{x)T^PR{y)aRQ{z) 
+GGPQR{Pa,Pb) aQp{x)apR{y)t^RQ{z) 



+G7PQR{pa,Pb)crQpix)apR{y)aRQ{z) J . (3.11) 

Tlie divergent parts of tfie one-loop functional (jH.lip are contained in the functions A, B{0) 
defined in Eqs. ()B.1|B.14|1 . Renormalization amounts to absorbing those divergences in the low- 
energy constants Lj of the next-to- leading Lagrangian ()C.1|) . rendering those constants scale 
dependent at the same time. In other words, to obtain the renormalized, scale-independent 
generating functional Z4 in Eq. ()3.1|) . the divergences A(/i) defined in ()B.12|1 are to be dropped 
in Z^^^, with the Li being replaced by the renormalized constants L^{fi) in the local func- 
tional Z^'^^. The analogous procedure will be understood for the nonleptonic weak part of the 
generating functional of discussed in the following section. 

The momentum flow shown in Fig. ^was chosen to fulfill the following criteria: 

• The generating functional should contain as few terms as possible. 

• The numerator of the three-point functions, i.e. (/c^ — M^)((/c— pa)^ — M^)((A;— Pf,)^ — M^), 
should be symmetric under the interchange of Pa and pf,. 
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4 Nonleptonic weak interactions 



To describe \AS\ = 1 nonleptonic weak processes, one has to insert the nonleptonic weak 
Hamiltonian in the appropriate Green functions. At the hadronic level, this amounts to includ- 
ing exactly one vertex from the effective nonleptonic weak Lagrangian in the various parts of 
the generating functional. 

At leading order, the nonleptonic weak Lagrangian of 0{Gfp'^) describing \^S\ = 1 pro- 
cesses is given by [13] 



^2 



+ h.c., (4.1) 



where L^ij denotes the zj-component of the matrix L^. The weak mass term proportional to 
{X{Wx + X^U)) can be transformed away by a suitable field redefinition [6, 14]. The octet and 
27-plet coupling constants Gs and G27 can be extracted from K ^ nn decay rates [15]. We 
work consistently to first order in Gs and G27- 

To 0{p^), with terms linear in Gs,G27 included, the generating functional is still given by 
the classical action ()2.1U|) with 

C2 = Ci + CY . (4.2) 

The classical solution U now satisfies a modified EOM where the weak part due to is added 
to piTl) . 

At O(p^), the structure of the generating functional is again given by Eq. 1)3.11) . The tree- 
level part Z^'^^ now contains the action for the nonleptonic weak Lagrangian of 0{GfP^) 
reproduced in Eqs. ()C.4IC.6|) . The chiral anomaly affects weak amplitudes both through the 
weak part in the classical solution U in Z^^isn and intrinsically through terms proportional to 
the e-tensor in the weak Lagrangian in ()(I3jl [16]. 

It remains to determine the weak part of the one- loop functional Zf^^. The operator 5 
defined in Eq. fl3.5|) is now of the form 

5 = 5s + 5w . (4.3) 
The weak perturbation 5\y is given to first order in Gp by the following expression [6]: 

5w = {f J^, + {f J^, n + , = -1 at; + 1[T - Tt, f,] , (4.4) 

a^ = u + ^[d^, N^] - i[rf^ [d^, a]] - a} + Ta + aT^ , 
Ta, = ^ffu^ ' d,A = d,A+[t^,A], a = a-a\ 
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The matrices (T° and are obtained from a and a by switching off the external fields. In the 
weak perturbation there are no terms without external fields because a diagonalizing transfor- 
mation similar to the one in Ref. [17] has been performed. The components of the antisymmetric 



matrix A^,, are 



1 1 



Nf,ab = -^^^{{K32,y^^}[Xa,Xb]) - -K{Ks2iXayf,Xb- Xbyi^Xa)) (4.5) 

-C{[\a,\b]K){Ky^) + ^C{K[K,y,]){K\b) - ^C{MK,y,]) {KX^) + h.c. 
The coefficient k, stands for 

K = AG8F'^ (4.6) 
and the matrix K is defined by matrix elements 

Kij = uXijv) with {\ij)ki = SikSji . (4.7) 

The operators ({KO){KV) are due to the 27-plet contribution and are defined as 

C{KO){KV) = AG27FH,,,u{K,,0){K,iV) , (4.8) 

with coefficients 

_ 1 _ _ 1 

^23, 11 — ^11, 23 — 2 ' ^21,13 — ^13, 21 — - , 

tij,ki = otherwise . (4.9) 
The expressions for the remaining symmetric matrices are 

+ ^-C{Xa[K,y,]){KXb) + ^C{Xb[K,y^]){KXa)+h.c. , (4.10) 

1 1 

(^ab = --K,{{Xa, Ab}{y^ -ft'32}) + -^i^{K32 {Xay'^Xb + Xby'^X^) 

+^i^{{K2.2,y''}{XaypXb + XbypXa)) - ^t^^y'' K2,2y p{Xa, Xb}) 

-^C([i^,/]A„)([ir,l/,]A,) - \aKyP){{K,yp}{Xa,Xb}) 

+C{K{Xay'Xb + Xty'Xa)){Kyp) +h.c., (4.11) 

aab = \'^{{Xa,Xb}Ks2) + ^C{XaK){XbK) +h.c. (4.12) 

To get the weak functional one must replace in Eq. (jH.llj) f ^ by f and a by once in each 
term and insert the appropriate binomial factors. 
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5 Mathematica® program Ampcalculator 

The program Ampcalculator [18] evaluates the complete amplitudes for strong and nonleptonic 
weak mesonic processes to O^p"^) and to 0(0^) (counting an external photon as 0(0^)), with 
the following exceptions. 

• In the Wess-Zumino-Witten functional [12] incorporating the chiral anomaly, only terms 
with explicit spin-1 fields are processed. In other words, we omit for practical reasons^ 
the part that cannot be written as the integral of a finite polynomial in U, d^U in four 
dimensions but can be written as a five- dimensional integral, with Minkowski space as 
boundary of the five-dimensional manifold. 

• For radiative semileptonic decays (e.g., for l^ui'-f), the Bremsstrahlung ampli- 
tude(s) for the photon(s) coupling to the charged lepton must be added by hand. For 
instance, for a single photon in the final state the program only calculates hadronic matrix 
elements of the type 

M^, = i j (i^xe^^-^(hadrons|ry;'°^(a;) (K(0) - A,(0))l^^l=' |M(eson)) (5.1) 

and then multiplies the result with the appropriate factors to get the (partial) amplitude 
for the decay M — hadrons + / + i/; + 7. 

By default, the program reduces the one-loop amplitudes down to the basic loop functions 
A,B,C defined in App. |B] In certain cases, the result may be more transparent and even more 
suitable for numerical treatment without performing the recursion relations for the various 
loop integrals collected in App. |Bl The user has the option to skip the recursion relations by 
setting the variable norecurs=l. 

For running Ampcalculator, a list of external mesons, photons and must be specified 
together with their momenta in an input file (cf. App. All particles are assumed to be 
incoming. The program checks for charge conservation and exits if the total charge of the 
listed particles is non-zero. For convenience and better readability of the output, the user 
can also provide two lists with scalar products of momenta and polarization vectors (also 
used for virtual photons for simplicity). The first list contains the scalar products themselves 
and the second one the corresponding scalar variables. In this way, momentum conservation 
can be implemented. The lists can also be used to redefine momenta (see Sec. 15.11 for an 
example). For semileptonic transitions, the program calculates the V — A hadronic amplitude 
and then multiplies with the appropriate factors to get the full amplitude (except for lepton 
Bremsstrahlung) including the lepton matrix elements or Z^, where 

= u{pi-)-ff,{l--^5)v{p^^) . (5.2) 
^The associated vertices contain at least five meson fields. 
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The masses of the pseudoscalar mesons (always in the isospin hmit) need not be specified 
exphcitly for the squares of meson momenta. However, the user should specify (for his/her own 
convenience) if the external photon(s) is (are) real: for an on-shell photon with momentum k 
and polarization vector e[k], the lists should contain the declarations for k ■ k = k ■ e[k] = 
(cf., e.g., Sec. 15. 3p . In the amphtudes of 0{p'^), by default is replaced using the Gell-Mann- 
Okubo relation 

M', = I (4M| - M'^) , (5.3) 

except as argument of a loop function where for better readability is always kept. The user 
has the option to keep everywhere. The complex conjugate of a quantity G (e.g., Gg, iVj, 

Vus, • • • ) is denoted generically as G. 

Running Ampcalculator may be quite time consuming, especially for nonleptonic weak 
processes, the time increasing strongly with the number of external states. However, the big 
advantage is that it is computer time rather than the physicist's time that is being consumed. 
Especially for exploratory purposes, the program provides options for calculating the tree 
amplitudes (at lowest order or to 0(p^)) and the loop amplitudes separately. Details can be 
found in App.^where the input file for Ampcalculator is listed. The input file ampcalculator. nb 
and the subroutine amp calculator sub. nb are available for general use [18]. The results are 
written into a separate output file. 

To perform CHPT calculations, another Mathematica® package, called PHI [19], exists. 
That package makes use of the more general package FeynCalc [20]. In contrast to Ampcalcu- 
lator, FeynCalc/PHI is not a fully automatic program. In order to use FeynCalc/PHI, detailed 
knowledge of Mathematica® and of the packages FeynCalc and PHI is needed. Another dif- 
ference between FeynCalc/PHI and Ampcalculator is the way the amplitudes are calculated. 
FeynCalc/PHI generates and uses Feynman rules whereas Ampcalculator performs a functional 
differentiation of the generating functional. Both packages are completely independent allowing 
for an excellent check of the results. 

In the remainder of this section we present input and output for three examples. Although 
the examples are nontrivial from the point-of-view of computer time needed the amplitudes 
are relatively simple and can easily be compared with results available in the literature. The 
indicated execution times were obtained running Mathematica®5 on a PC with 512 MB RAM 
at 1.6 GHz. 

In addition to the three examples shown below, we have checked the existing results for the 
processes 77 Tr+Tr" [21], e+e" An [7,22], K+n+ K+Tr+ [23] and 377° [24,25]. 

5.1 The semileptonic decay 7r^/+z// (Kj^) 

The program Ampcalculator evaluates the hadronic V — A amplitude and then multiplies it 
with the appropriate factor to get the full decay amplitude. With the external momenta defined 
through 
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with 



q = Pi+Pi 



p- p', p = p + p' , 



(5.4) 



a possible input for the hsts ExternalParticles, replistl and rephst2 is given below. Of course, 
the choice of momenta and kinematic variables is up to the user. Remember that the program 
interprets all particles as incoming. 



ExternalParticles = {K+[p], W_[Q], 7ro[r]}; 

replistl = {sp[Q, Q], sp[r, Q], sp[p, Q], sp[p, r], p^^, r^', Q'^}; 

replist2 = {t, (M2,-M2-t)/2, {-Ml+Ml-{)/2, (-M|,-M2+t)/2, (P'^+q^)/2, (-P^+q'^)/2, 

The complete result of O(p^) is printed out in a separate output file (execution time: 
0.30 hours). 

Complete amplitude to 0(p''4) 



(t - 5 Mk^ - Kjr^) (P^ t + q'* (-Mk^ + M^2)-) 
192 7r2tF^2 

(_2q/^ (Mk^-M^2) (t-4MK2 + 2M^2-)^p/..,. (t - SMk^ + 4M^2)) AfMs^] 

8tF,2 (M,2_M,2) 

(P^ t (t - 4Mk') + 4 q^ (t + Mk') (Mk' - M^^)) AfM^^] 
16tF^2 (Mk2-M,2) 

3 (P^t (t - 4Mk2) + 4q^MK2 (Mk^ - M^^)) ^[M^^] 

letFT^MT^MT) 



Gf1/.V, 




P'' 2 q'^ (Mk^ - M^2) 



+ 



Lr (pMt + q^ (-Mk' + M/)) 
F 2 




144t2F^2 



The result agrees with Ref. [26]. 



13 



5.2 The nonlep tonic decay K — ^ tt tt^ 

In this decay, the complex conjugate couphng constants G27 and i?[ appear in the amphtude. 
The assignment of momenta and scalar products is straightforward. 



External Particles = {K_[pi], 7r+[qi], 7ro[q2]}; 



replisti = {sp[pi, qi], sp[pi, qs], sp[qi, qa]}; 



replist2 = {-i Ml,, -i Ml,, -M^+iMI,}; 
Complete amplitude to 0(p''4) 



iG27l -F^ (Mk'-M/ 



80 (2Mk^-Mk'M^^ -M/) 
3F^ 



20 (Mk^ + 2 Mk^ _ 3 M^4-) ^ 5 ^j^^4 _ 3 ^^2 ^^2 ^ 2 M^^) 



3F. 



96 7r2 



3F, 



3 F^ 



10 (Mk^ - M,2) (2Mk2 ^ (Mk' - K^) K2 



3 Ftt 3 F^ 

5 (3 Mk^ + M^2) a[Mk2] 5 (4 Mk^ - 22 Mk^ M^^ + 29 M^^) AfM^^j 



12 F, 24F,M,2 
5 A[M^2] 5 ^^^4 _ 3 ^^2 ^2 ^ 2 M^^) B[Mk2, M^^ M^^] 



8 F^ _ 6 F^ 

5Mk^ {-n^ + K^) B[M^^Mk2,M^2] 



72 

5Mk" (5Mk^- 13Mk^M^^ + 8M/) B[M^^ Mk^] \ 
24F^M^2 J 

With the appropriate renaming of coupling constants, the result (execution time: 1.48 
hours) agrees with previous calculations in Ref. [25,27]. 



5.3 The nonleptonic decay — ^ tt 77 

Finally, we present the amplitude for a decay where three-propagator loops contribute: 

K'ip) ^ 7r°(p')7(?i)7(?2) . 

In the limit of CP conservation, the even-intrinsic-parity amplitude determines the amplitude 
for Kl — > 7r°77 whereas the part with the e-tensor yields the corresponding Ks amplitude. 
The following assignments (all particles incoming) correspond topi = —p', ki = —qi, k2 = — ?2- 
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External Particles = {Ko[p], 7ro[pi], -fih], -f[k2\}; 

replisti = {sp[p,pi], sp[ki,k2], sp[ki,e[ki]], sp[ki,e[k2]], sp[k2,e[ki]], sp[k2,e[k2]], sp[e[ki],e[k2]], 
sp[p,e[ki]], sp[p,e[k2]], sp[pi,e[ki]], sp[pi,e[k2]], sp[ki,ki], sp[k2,k2j}; 

replist2 = {| (s-M^^-M^), |, 0, -qr 62, -q2- ei, 0, ey €2, p- ei, p- 62, -p- ei+qa- ei, -p- e2+qr 62, 0, 0}; 

Complete amplitude to 0(p''4) 

f (2 qi • €2 q2 • ei - s ei • £2) Mk^ 

+-V^C[0,0,-^,M^^M^^M^'] (2qi-e2q2-ei-sei-e2) M^' (s - M^') 

+ fM2^ M 2N ^C[0, 0, MK^ MK^ Mk'] ("2 qi • £2 q2 • 61 + s 61 • 62) Mk' 



X (sMk'-M/-sM^' + M/)) 
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(2 qi • 62 qs • ei - s ei • 62) (-13Mk^ + 15 (s - M^")) 



+ 2n/2C[0, 0, M,^ M,^ M,2] (2 qi • 62 q2 • ei - s • £2) M,^ 
X (-5M/ + 24M^2 (-S + M^') + Mk' (q s + M^')) 

+ I m2n 2\/2 C[0, 0, MK^ MK^ Mk'] (-2 q, . 62 q2 • + s • 62) Mk' 
o s [i\ — iv j z 

X (21 Mk^ + 6 M^' (s - M^2) + Mk^ (-21 s + 5 M^' 

1. , pV2 6^^'--G8 (s-Mk^) (Mk^-M^^) (pgki,-kigp,J (6(kO)^(6(k2)), 

2v^e«'"^-G2T (s-Mk^) (Mk2-M^2) (p^ - k^^ p. ) (e(kO), (e(k2)), 




7r2 (s - M^2) (3 s - 4 Mk^ + M^^) 

The amphtude (execution time: 1.68 hours) agrees with what is available in the literature 
[28,29]. The function F defined in Ref. [28] is related to the three- propagator loop function as 

Fis/M"^) = 1 + 327r2 C(0, 0, -s/2, M^, M^, M^) . (5.5) 

In the anomaly contribution, was replaced via the Gell-Mann-Okubo relation. To our 
knowledge (and for good reasons), the 27-plet kaon loop amplitude and the anomalous contri- 
bution proportional to G27 have not been published before. 
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6 Conclusions 



The generating functional is a convenient quantity for a compact representation of Green 
functions and S-matrix elements. Especially for an effective field theory like chiral perturbation 
theory, with its many coupling constants and derivative couplings, the standard diagrammatic 
calculations can be quite cumbersome and they must be performed for each process separately. 
The great advantage of the generating functional is that the renormalized amplitudes can be 
obtained once and for all. 

Extending previous work, we have presented in this paper the generating functional of 
Green functions for chiral SU(3) for both strong and nonleptonic weak interactions in the 
meson sector to next-to-leading order in the low-energy expansion and for at most six external 
states, a photon counting as 0{(fP'). Such processes require the inclusion of one-loop diagrams 
with up to three meson propagators and they comprise almost all strong and nonleptonic weak 
mesonic transitions of physical interest. 

The main purpose of this work is to reproduce and to check previous calculations in a 
simple and straightforward way but the functional presented here has also been used to derive 
new results [10]. Although the representation is compact the actual extraction of a specific 
amplitude via three-dimensional flavour traces can be quite laborious, the toil increasing of 
course with the number of external states and especially for nonleptonic weak transitions. 
We have therefore provided as an essential part of this work the Mathematica® program 
Ampcalculator for general use that performs the necessary manipulations upon input of the 
relevant external particles with their momenta. 

We have described the necessary input for the calculation of three semileptonic and nonlep- 
tonic K decay amplitudes and presented the resulting output. As far as available, the results 
agree with previous calculations. 

The generating functional presented here can be extended in several ways. In addition 
to radiative corrections for strong processes [8,9], photons and leptons can be included as 
dynamical degrees of freedom for similar treatments of radiative corrections for semileptonic 
[30] and nonleptonic weak decays. 
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A Input for the program Ampcalculator 



We list here the contents of the input file ampcalculator. nb. During execution of the input file, 
the subroutine ampcalculatorsub.nb is called where the actual computations are performed. 
The result is then printed out in a separate output file. 

Program Ampcalculator 

The program generates the amplitude for any strong or nonleptonic 
weak process of pseudoscalar mesons, external photons and W bosons to 
next -to-leading order in the low-energy expansion (chiral SU(3)) with at 
most six external particles (a photon counts as two particles) , at most one 
W (semileptonic decays) and with at most three propagators in the loop. The W 
is assumed to couple to a lepton pair indicated in the output by the leptonic 
matrix elements 1^ and 1^, defined by 

V = ^(p..)7m(1 - 75)v(pi+) , 1^ = u(pi-)7/,(l - 75)v(p^) , 

where , pi are the momenta of the neutrino and of the charged lepton. For 
radiative semileptonic decays, Bremsstrahlung off the charged lepton is not 
included and must be added by hand. In the anomaly functional, only terms 
containing external gauge fields are included. 

Select up to six particles (a photon counts as two particles) of the list 

7r_ , TTo, , , Kq, Kq, rjs, 7, W+, W-} . Assign a momentum to each 
particle. In the list ExternalP articles all particles are supposed to be 
incoming. Define the scalar products sp[pj, pj] that you want to replace by 
polynomials of kinematic variables in the lists replistl and replist2. One should 
not use the capital letters F, G, L, N and R in the lists replistl, replist2 
and ExternalP articles as they are reserved for coupling constants. Small 1 is 
reserved for the leptonic matrix element. For semileptonic decays the index 
should only be used for momenta that are contracted with the leptonic matrix 
elements 1^ or 1^. Isospin conservation is assumed. 

The program terminates if either 

i) some of the external particles are not contained in the list above, 
ii) charge is not conserved, 
iii) there are more than one external W, 
iv) strangeness change |AS|>1. 

Notation: 

A[M^] is the renormalized single-propagator loop integral 
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A[M^] = — - — — log(M^//i^) , where M is a meson mass and is the renormalization 

(47r) 

scale. B[s, M^,M|] is the standard two-propagator loop function subtracted at 
s=0, i.e. B[0, , M|]=0 . All loop functions and recursion relations can be 
found in App. B. Gg, G27 are the nonleptonic weak LECs of order . L^,N^,R^ 
are the strong, weak octet and 27-plet LECs of O(p^), respectively, all 
renormalized at scale ji. The corresponding Lagrangians are reproduced in App. 
C. Gs, G27, Ni, Ri, are complex conjugate LECs. 

Options : 

Set the following variables equal to 1 to get partial results 
(otherwise the variables should be set to 0) : 

onlytreep2=l amplitude of O(p^) 

onlytreep2and4=l tree level amplitude up to O(p^) 

onlyloops=l one-loop amplitude only 

By default (nogmo=0) , the Gell-Mann-Okubo relation is applied in the 
amplitudes of O(p^) to express in terms of M]^, . For better 
readability, is never replaced as argument in loop functions. Setting 
nogmo=l, is kept everywhere. 

If the recursion relations for the loop integrals should not be performed 
set norecurs=l. 

onlytreep2 = 0; 
onlytreep2and4 = 0; 
onlyloops = 0; 
nogmo = 0; 
norecurs=0; 

(* Example: *) 
ExternalParticles = {ttq [pi] ,7ro [^2] ,7r+ [ps] ,7r_ [^4] } ; 

(* scattering amplitude for ttq [pi] + ttq [^2] — tt- [— P3] + 7r+ [— P4] *) 
replist 1 ={sp [pi ,p2] , sp [pi ,p3] , sp [pi ,p4] , sp [p2 ,P3] , sp [p2 ,P4] , sp [ps ,p4] } ; 

replist2 ={i(s-2M^), ^{t-2M',), ^ - 2M^) , i - 2M^) , ^{t-2M',), 
(* *) 

(* Insert the correct path in the next line.*) 
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notesub=NotebookOpen ["C : \Math5\prograins\Unterdorf er\ampcalculatorsub . nb"] ; 
SelectionMove [notesub , All , Notebook] 



SelectionEvaluate [Notebook] 



B Loop integrals and constituent functions 

The standard functions A, B, B^, C and Cij occurring in one-loop integrals with up to three 
propagators are defined through the following relations [Cp denotes the Feynman integration 
contour) : 

d<^k 1 



I Jc 



Cp {2t:Y (A;2 - Ml) ' 



(B.l 



^ ^ iJc^{27iY {k^~Ml){{k-p)^-M^)' ^' 

Bip') = {1} , (B.3 

{k,}^p,Bn{p') , (B.4 

{kfxkS\ ^ gi,uB2o{p^) +P^,PuB22{p'^) , (B.5 

nx\\ = -[ - fB6 

i Jc, {2nY {k^ - Ml){{k-paf - Ml){{k-p,f -Ml)' ^ ' 

C = {{1}}, (B.7 

= Pafi Cua + Pbfji Cllfe , (B.8 

{{k,j.ky}} = Qi^u C20 + PatxPai^ C22a + PbuPbu C226 + {paixPbu + Pbu,Pau) C22 , (B.9 



{ { k^ k^ kx}} = Paf^PauPaX CsSa + PbfiPbuPbX C?,?,b + 

+ {Paf,PauPbX + Pa^iPbuPaX + PbpiPauPax) C^Sa + 
+ {PbfiPbuPaX + PbfiPauPbX + Pai^PbuPbx) C'sSfo + 
PauQ^X ^" Pa 

+ {.Pbn9uX + PbyQiJiX + PbxQuv) C316 , (B.IO) 



with = Cij{pl,plpa ■ Pb, Ml, M2 , Ml), Qjb = Cij{plpl,pa ■ Pb, Ml, M^ M2) and = 
Cij a = Cij I,. The functions A and B can be split into a finite and a divergent part for d ^ 4: 

A{Ml) = A{Ml) - 2MlA{ii) , Bip^) = Bip^) + B{0) , (B.ll) 
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where 

^^'^^) = -(4#^"#' ^^^) = ^(l3I-2(M47r) + r'(l) + l)), (B.12) 
^^^^ ^ 32;^i^+ p2 ^^Mf "M|-M^^"m| 

2 \2' 



(p2 + ^A(p2,M2,M2))2-(M2 -M^ 

X m 



and 



A(x, y, z) ^x'^ + + z'^ - 2{xy + yz + zx) , (B.13) 



The coefficient functions are of the form 



i?20(p') = 



2p2 



288 7r2 6 



2N F 



288 ttV 3p2 



(B.15) 



The pure C function is given by 

C(ApI p„ . p. Ml. Ml MI) = -^-i. g £ ^L, ^-^ ) - L, 



1 1 ^ ^ ^ ^ ( -r • ( 7 \ -r • ( 9 1 



(B.16) 



with the dilogarithm 

U2{z) = - -ln(l-t) (B.17) 
Jo t 

and 
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V 



1- o ^ ± 

Pa 



pI 



PD=Pa-2Pa-P6+P6, KpI^pI^Pd) 



(B.18) 



The quantities X2, -zf and X3, are obtained by cyclic interchanges of {pi^VDiPD ^iid 
(M|,M|,M|). 

In the following 3) is defined as M|,, M|) and 5(2, 3) as 

B{p\ — 2pa ■ Ph + pI-i Mq, M^). The same notation is valid for Bu, B20 and i?22- 



C'lla — 



HllaPb^ - HiibPa-Pb 
pJPb^-Pa-Ph'^ 



1 , 2CM| + E(2,3)-Cna/i-giu(Mp^-M/+p,^ 



C20 — -^-j — J + 

64 TT^ 

H21aPb' - H22bPa-Pb 



C22a 



Pa^ Pb^ - Pa ■ Pb^ 



G 



'22 



H22bPa^ - H21aPa-pb 
Pa^ Pb^ - Pa ■ Pb^ 



_ HsiaPb"^ - H32bPa- Pb ^ _ H32bPa^ - H^iaPa ' Pb 



Pa^Pb^ -Pa-Pb^ 



Pa^Pb^ -Pa-Pb^ 



^ _ HsOaPb^ — HsObPa ' Pb 
^31a - 2^2 



Pa^Pb-Pa-Pb 



2 ' 



_ -i?(l,3)+i?(2,3) + C/i _ i?(2,3) + Cna/i-i?ii(2,3) 

-nila — 7, 1 -"21a — 7, ^-^20 ; 



„ _ Ciu/i-Sn(l,3) + 5n(2,3) ^ _ C20 /i - ^20(1, 3) + ^20(2, 3) 

-"22a — 7 , -"SOa — 



H: 



2 ,0.^ 2 

5(2. 3) + h C-na 2 Z?n(2. 3) + Z?.2(2. 3) 



31a 



— 2C31a , 



32 a 



/i C'226 — -622(1, 3) + -622(2, 3) 



(B.19) 



The results for the integrals agree with those in Ref. [31]. 



21 



The constituent functions of the one-loop functional are of the following form: 

GtpQip) = P'P" M'p, Ml) + 5^^^ h{p\ Ml, Ml) , 

G^2PQip) = ^P'' {^B{p\Ml,Ml)-B,,{p\Ml,Ml)^ , 
G^pq{p) = -B{p\mI,MI), 



GTpQRiPa,Pb) = i (caPa^Pa.Pa^ + ChPh^PKPb^ + daPa^Pa.Pb;, 
+eaPa^Pb,Pa^ + faPb^Pa^Pa^ + dbPb^Pb,.Pa^ + fbPb^Pa.Pbx 
+ebPa^Pb,Pb^ + QaPay^ + haPa^Q^^ + iaPa^Q^'^ 

+hbPby^ + gbPb.g"^ + ibPby) , 

G5PQR{Pa,Pb) = Jag^'" + KPa^Pa^ + hPb^Pb^ 
+LPa^Pb,.+maPb^Pa. , 

GQPQR{Pa,Pb) = i {naPa^ + UbPh^) , G7PQR{pa,Pb) = \ g , (B.20) 

with for example 

Ca = c{pl,plpa ■ Pb, Ml, Ml, Ml) and Cb = c{plpl,Pa ■ Pb, Ml, Ml, M^) , (B.21) 
where 



a{p\Ml,Ml) 



4(^n(p^ Ml, Ml) - ^22(P^ Ml, Ml)) - B{p\ Ml, Ml) 

4 



h{p\ Ml, Ml) = -B2o{p\ Ml, Ml) , 

_ C*!! a ~ 4 C22 a + 4 C33 „ _ Ciia — '2 C22 a ~ 2 C22 + 4 C33 „ 

_ -C/2 + 2Cna + Cllb - 2C22a - 4(722 + 4(733,, 

_ -2 (C'22-2C'33a) _ -2 (C20-2C31a) , _ 4 C31 g 

/a 3 -I ga 2 ) "-a 3 ) 

. _ -2 (C20 - 2C3ia) ^ _ c.(^ h - r or 
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c 

la — + Cua + C'llb — 2 C22 



m„ 



-2C. 
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^ — 2 Ciia 



(B.22) 



C Lagrangians of O(p^) 

The strong Lagrangian of O(p^) is given by [3] 

£f = L^iD^U^Df^Uf + L2{D^U^D,U){D''U^D''U) 

+L7(x^f/ - xf/^)' + U{x^Ux^U + xf/^xf/^) 

-iL,{F^'D^UD,U^ + F^D^U^DM) + L^o{U^ F^^U Fl^,) (C.l) 

with 

-^i = ^^^^i^ - ^^^^/^ - i [^M' ^i^] ' -^i = ^^^^^^ - ^^^^^^ - i [li^^ lA ■ (C.2) 
The nonleptonic weak Lagrangian of 0{GfP'^) consists of an octet and a 27-plet part: 

CY = Cf' + Cf^' . (C.3) 

Neglecting terms that would only contribute at 0{G'^) (with external fields in addition to 
the nonleptonic weak transition), the = 1 octet Lagrangian can be written as [6,32] 

+N5 {X{S,Df,U^D^'U}) + N6 {XD^Jj'^U){U'^ D^'U S) 
+Nr {XS){D^U^D''U)+Ns {XD^U^D''U){S) 
+N, i{X[P, D^U^Di^U]) + N^o {XS'') 

+Nn {XS) {S) - N,2 (XP^) - A^i3 ( AP) (P) (C.4) 

+A^i4 ^(A{y^^ D^U^'DM}) + A^i5 i{XD^U^UV''''U'' DM) 

-7Vi6 ^(A{A'^^ D^C/t - Nir iiXD^^UA^^M^DM) 

+2iVi8 {XiF^U^FR^M + U^Fr^MFD) 

+N28 ie^,.p^{XD^U^U){U^ DMDM^ DM) 

+2iV29 e^,p„{X[U^Fi;,MM'U^DM]) 

+iV3o e^..p„{XU^DM){V''''DM^U) 

-N^i e^,upa{>^U^DM){A'"'DM^U)} + h.c. 
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with 



S = x^U + U^x, P = i{x^U-U^x), 
V,. = U^F^^U + f;:,, A,, = U^F^,U-F^,. (C.5) 

For the 27-plet Lagrangian we use the form given in Ref . [33] : 

£^^,27 ^ G27FHij^klIji,lk + h.c. , (C.6) 

+R,S,,{L^L,),, + R,{S,L,}^^{L^),, 

+i?9i[p, L,\^m,i +Rio{s) {L,\^ m,, 

+Rn S,, Ski + R12 P^J Pki + i?i3 i (f^-),, {L^ L''),^ (C.7) 
+Rui {a,.),^{L^n,, + Ri,i K., L^\^ in^i 

+i?i6i K., L% (1^,1 +R17 {vn, 



•-HU-, ^ \ij V-^ )kl ^ -"^17 \^iiv}ij \^ )kl 

+Ris M,, {a'^n,i + Ri9 (a^.),, {a^ki 

+R20 i e^'^^" {L^ L,)^^ {Lp L.),^ + R21 i e^^^'^ (L^)^^. (L, L, L,)^^ 

+R22 6^^^- {L,)^^ (L. v,^),, + R2, e^'^P'^ (L,)^^. (L. a,.),, . 



The coefficients iij^ki are defined as 

_ ~ _ ~ _ ~ _ 1 

il2,31 — ^31, 12 — ^32,11 — ^11,32 — g , 

_ 1 

^22,32 — ^32,22 — ^32,33 — ^33,32 — "77 ) 

6 

iij,ki = otherwise . (C.8) 

The low-energy constants Lj, Ni and i?j are dimensionless. Moreover, the Lj are real whereas 
Gg, G27, Ni, Ri can be complex in the presence of CP violation. 
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